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INTRODUCTION 


Prominent  among  the  materials  used  in  modern  solid  state  technology 
are  the  III  - V and  II  - VI  semiconducting  compounds.1-78*  These  substances, 
typified  by  gaVlium  arsenide,  are  most  often  of  the  zinc-blende  structure, 
crystal  class  43m.  Much  work  has  been  done  on  the  growth  and  characteriza- 
tion of  these  crystals,  because  of  their  practical  importance,  to  the  point 
where  high  purity,  high  resistivity  materials  are  becoming  a reality.  These 
compounds  are  piezoelectric  in  nature,  but  the  effect  is  "short-circuited" 
if  the  resistivity  is  too  low.  The  piezoelectric  effect  couples  the  elec- 
tric and  acoustic  fields  present  within  the  crystal  and  makes  possible  de- 
vices that  depend  upon  the  interaction.1’2 

The  theory  of  acoustic  wave  propagation  in  piezoelectric  semiconducting 
crystals  is  due  to  Kyame 8 , and  to  Hutson  and  White4  . Arlt5  gave  an  ap- 
proximate analysis  of  semiconducting  vibrators  for  the  purpose  of  finding 
the  piezoelectric  coupling  factor.  Since  that  time  a number  of  advances 
have  been  made  in  the  theory  of  piezoelectric  insulators,  that  have  not  been 
extended  tc  cover  the  case  of  semiconductors6  . These  include  the  solution 
to  the  piezoelectrically  coupled  simple  thickness  modes  of  a plate  vibrator, 
the  use  of  doubly  rotated  cuts  for  resonators,  synthesis  of  exact  equivalent 
networks  using  transmission  lines,  and  others. 

In  this  report  we  undertake  to  explore  some  of  the  consequences  of_acou- 
stic  wave  propagation  in  piezoelectric  semiconductors  of  crystal  class  43m 
from  the  vantage  point  of  the  advances  mentioned  above.  Doubly  rotated 
plate  vibrators  are  considered,  as  are  coupled  modes  and  transmission  line 
equivalent  networks.  It  is  also  shown  how,  by  using  such  crystals  as  pie- 
zoelectric resonators,  sensitive  measurements  of  oxide  films  may  be  made  in 
a nondestructive  manner. 

By  way  of  providing  practical  examples  of  the  properties  of  doubly  ro- 
tated resonators  of  cubic  materials,  calculations  have  been  made  using  gal- 
lium arsenide.  This  material  is  widely  used  because  its  wide  bandgap  pro- 
vides for  high  temperature  operation  and  its  high  mobility  permits  high  fre- 
quency capability. 


PLANE  ACOUSTIC  WAVE  PROPAGATION  IN  CUBIC  CRYSTALS 


The  elastic  stiffnesses  yc.A/t  , piezoelectric  constants  elA  , and  di- 
electric permittivities  of  insulating  crystals  depend,  for  their  sym- 

metries, on  the  crystal  class.  For  crystals  in  cubic  classes  4”3m  and  23, 
the  elastopiezodielectric  matrix  is  given  in  Table  1.  It  is  seen  that  these 
materials  are  dielectrically  isotropic;  in  addition,  they  have  only  one  pie- 


* See  list  of  references  beainnina  on  Daae  dn 


zoelectric  constant,  eIH  , which  reduces  the  transformation  of  the  for 
rotated  axes  to  a problem  in  simple  cjeometry. 

For  coordinate  axes  X-'  transformed  by  an  orthogonal  matrix  <X(.j  from 
the  crystal lographic  axes  X;  we  have 


In  Table  2 the  left-hand  entries  are  the  subscripts  ( k/x)  representing  the 
quantities  ( e^/eIH  ),  the  normalized  transformed  piezoelectric  constants; 
the  numbers  in  parentheses  preceding  the  constants  are  numerical  factors. 

On  the  right-hand  side  are  products  consisting  of  ckjj  components;  only  the 
( ij  ) subscripts  are  written.  Table  2 follows  from  writing  out  the  trans- 
formation law  for  tensors  of  the  third  rank,  and  applying  the  simplifica- 
tions in  the  e;x  from  Table  1. 


When  the  transformation  of  Equation  (1)  is  expressed  in  terms  of  the 
angles  4>  and  & used  to  describe  doubly  rotated  cuts  of  orientation7 
( YX  url)4>l  9,  the  result  is  the  further  simplification  shown  in  Table  3. 
On  the  left-hand  side  are  the  indices  representing  the  transformed  quanti- 
ties ( eK>u/eIH  ),  and  on  the  right-hand  side  are  the  angular  dependencies 
of, these  quantities;  numbers  in  parentheses  denote  equalities,  e.g.,  €3^= 


Acoustic  wave  propagation  in  an  arbitrary  direction  in  crystals  may  be 
calculated  by  the  Christoffel  method6.  The  direction  of  plane  wave  propa- 
gation is  taken  to  be  the  X*  axis,  and  the  quantities  <XZ-L  in  Equation  (1) 
are  abbreviated  o( ■ in  the  following.  The  elastic,  piezoelectric  and  die- 
lectric constants  are  transformed  to  the  quantities  F1 , 3 , and  £ , re- 
spectively, by  the  relations  given  in  Table  4. 


When  an  acoustic  wave  travels  in  a piezoelectric  medium,  the  motions 
produce  internal  electric  fields  which  serve  to  increase  the  effective  elas- 
tic constant  governing  the  wave  motion.  This  effect  is  called  "piezoelec- 
tric stiffening";  the  stiffened  elastic  moduli  r are  given  by 


(2) 


The  P are  the  elements  of  a secular  determinant  that  determines  the  three 
stiffened  modal  elastic  constants for  the  three  types  of  plane  acoustic 
waves  that  can  propagate  in  the  assumed  direction  X'  . Associated  with 
each  mode  are  orthonormal  eigenvectors  i-  that  give  the  directions  of  par- 
ticle motion  associated  with  each  mode.  Tne  modes  are  ordered  according  to 


^ /cc  . (3) 

Mode  a is  the  quasi-longitudinal  mode,  with  particle  motion  nearly  along 
X'  ; mode  b is  the  faster  quasi-shear  mode,  and  mode  c is  the  slower  quasi- 
shear mode. 


The  semiconducting  nature  of  the  crystal  is  taken  into  account  by  per- 
mitting current  terms  other  than  displacement  and  piezoelectric  flows.  The 
current  density  due  to  drift  introduces  a conductivity  term3 


^ ej 
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TABLE  2.  ROTATED  PIEZOELECTRIC  CONSTANTS  FOR  CUBIC  CRYSTALS 

(1/6) • 1 1 = 12-13-11 
(1/6) -22  = 22-23-21 
(1/6) -33  = 32-33-31 


(1/2) - 12 

= 

22-23-11 

+ 

21-23-12 

+ 2 1-22-13 

(1/2) - 13 

= 

32-33-11 

+ 

31  -33-12 

+ 31-32-13 

(1/2) -21 

= 

12-13-21 

+ 

11-13-22 

+ 11-12-23 

(1/2)  -23 

= 

32-33-21 

+ 

31-33-22 

+ 31-32-23 

(1/2) -31 

= 

12-13-31 

+ 

11-13-32 

+ 11-12-33 

(1/2) -32 

= 

22-23-31 

+ 

21 -23-32 

+ 21-22-33 

14 

= 

11 "(22-33 

+ 

23-32)  + 

12- (23 -31  + 21-33) 

+ 

1 3 - (21  -32 

+ 

22-31) 

15 

= 

11 -(32-13 

+ 

33-12)  + 

12- (33 • 1 1 + 31-13) 

+ 

1 3 - (31  • 12 

+ 

32-11) 

16 

= 

11 -(12-23 

+ 

13-22)  + 

12- (13-21  + 11-23) 

+ 

13- (11 -22 

+ 

12-21) 

24 

= 

21 -(22-33 

+ 

23-32)  + 

22- (23-31  + 21 -33) 

+ 

23  • (21  -32 

+ 

22-31 ) 

25 

= 

21 -(32-13 

+ 

33-12)  + 

22- (33-11  + 31-13) 

+ 

23  - (31 -12 

+ 

32-11) 

26 

= 

21 -(12-23 

+ 

13-22)  + 

22- (13-21  + 11-23) 

+ 

23 - ( 11  - 22 

+ 

12-21) 

34 

= 

31  - (22  - 33 

+ 

23-32)  + 

32 - (23-31  + 21  -33) 

+ 

33- (21  - 32 

+ 

22-31 ) 

35 

= 

31 -(32-13 

+ 

33-12)  + 

32  - (33  - 1 1 + 31-13) 

+ 

33- (31  - 12 

+ 

32-11) 

36 

= 

31 -(12-23 

+ 

13-22)  + 

32 ■ ( 1 3 • 21  + 11-23) 

+ 

33- (11 -22 

+ 

12-21) 
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TABLE  3.  PIEZOELECTRIC  CONSTANTS  FOR  DOUBLY  ROTATED  CUTS  Cc 
CUBIC  CRYSTALS. 


11 

= 

0 

22 

= 

-3-sin 

20-cos20-sin  9 

33 

= 

-3-sin 

20 -cos  9-sii  "9 

12 

= 

+2-cos 

20 -cos  9- sin  9 

13 

= 

-(12) 

21 

= 

sin 

20-sin  9 

23 

= 

sin 

20-(2-cos20  - sin29)-sin  9 

31 

= 

sin 

20-cos  9 

32 

= 

sin 

20-(2-sin29  - cos20)-cos  0 

14 

= 

cos 

20-(cos^9  - sin^9) 

15 

= 

(31) 

16 

= 

(21) 

24 

= 

(32) 

25 

= 

(14) 

26 

= 

(12) 

34 

= 

(23) 

35 

= 

- (12) 

36 

= 

(14) 
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TABLE  4.  CHRISTOFFEL  MATRIX  FOR  CUBIC  CRYSTALS. 
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and  f for  classes  23,  432,  m3,  43m,  and  m3m:  5 for 

i < i 


classes  23  and  43m;  all  ■"  are  zero  for  classes  432,  m3, 

i 

and  m3m. 


while  a term  depending  on  charge  carrier  gradients  will  contribute  a diffu- 
sion term4.  If <r  is  finite,  the  effect  is  to  make  £ complex;  the  addition 
of  a diffusion  term  does  not  change  this  situation,  but  makes  the  relation 
more  involved.  Because  of  the  piezoelectric  stiffening  expressed  in  Equa- 
tion (2),  when  e becomes  complex,  so  does  cm,  and  the  modal  velocities 


where  P is  the  mass  density,  also  become  complex.  This  means  that  the  waves 
are  attenuated  as  they  progress.  The  piezoelectric  constants  are  only  in- 
fluenced to  a minor  extent  because  the  eigenvectors  become  complex,  and 
the  effective  piezoelectric  constant  for  mode  m is 

e~  = *-j  ■ 161 

If,  in  addition  to  modifications  toe  , viscous  forces  are  included, 3>6>8_ 
the  effect  is  to  add  a further  imaginary  term  to  the  elastic  stiffnesses  cm. 

The  effective  piezoelectric  constant  e^,  in  Equation  (6)  holds  for  elec- 
tric fields  directed  along  X*  . For  generation  of  surface  acoustic  waves9 
(SAW)  or  shallow  bulk  acoustic  waves1®  (SBAW)  excited  via  interdigitated 
electrodes9  , the  piezoelectrical ly  induced  mechanical  force  densities  Fj 
are  determined  from  the  relation 


(7) 


where  E*,;  stands  for  3EK/ay^  , and  the  uncontracted  indices  on  e have 
been  used.11  Table  5 gives  the  force  density  matrix  for  crystal  classes  h3w 
and  23,  using  the  e'  entries  from  Table  3 for  doubly  rotated  cuts;  only  the 
subscripts  are  shown.  One  may  see  that  the  synmetry  present  permits  the 
substitution  E w ~ E„  as  well  as  a number  of  other  relations. 

J'k  CR'J 


As  far  as  the  purely  acoustic  wave  velocities  of  SAW  propagation  on 
cubic  substrates  is  concerned,  explicit  formulas  can  only  be  given  for  pro- 
pagation along  the  cube  axis  and  along  the  face  diagonal  Ei-i °]  dir- 

ections. The  expressions  are:1*1 


For  propagation  along  [loo]  , with  the  abbrevation 

R = f or*  / /C„  , 


(8) 


( ± - yC ,,  R,  ) 

• ( l - R - /t,j  /yc,,1  ( ±~  R)  . 

(9) 

For  propagat 

ion  along  [no]  , with  the  further  abbreviation 

*.  = ( /C„  + /Cl2  +2  yc^)  h , 

(10) 

(l-  'Cjj  R 

( ( >€„  yC.  - >Ca  ) ~ (/Co//C„  - R)  . 

(ID 
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In  Equations  (9)  and  (11),  the  single  real  root  for  R yields  the  SAW  velo- 
city v from  Equation  (8). 


TABLE  5.  FORCE  DENSITY  MATRIX. 


El 

e2 

^3 

Ei , 1 

Ei, 2 

Ei  ,3 

E2.1 

E2 , 2 

E2 » 3 

E 3 » 1 
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E3 » 3 
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21 

12 

14 
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- 
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21 
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14 

12 

22 

32 

14 

32 

23 

-f3 

31 

14 

-12 

14 

32 

23 

-12 

23 

33 

Crystal  classes  43m  and  23. 


Table  6 presents  representative  values  of  elastic  constants  and  mass 
densities  for  selected  II  - VI  and  III  - V compounds,  along  with  the  anis- 
otropy ratios.  A,  defined  as 


A = Z / (/C,,  - /C(2 ) , (12) 

and  the  bulk  moduli,  defined  as 

8 - ( /C„  + Z /C  ,0/3  . (13) 


The  values  of  A and  B for  the  entries  in  Table  6 do  not  differ  appreciably. 
Because  of  its  high  technological  interest  gallium  arsenide  has  been  select- 
ed to  illustrate  some  of  the  foregoing.  Its  physical  properties  are  fairly 
representative  of  members  in  crystal  class  T3m.  Some  of  the  experimental 
values  characterizing  this  material  are  given  in  Table  7.  In  the  table, T*. 
is  the  first  order  temperature  coefficient  of  the  elastic  stiffness,  o<  <0 
is  the  first  order  thermal  expansion  constant,  and  the  /C>uV  are  the  third 
order  (nonlinear)  elastic  constants.  ' 

By  making  use  of  the  Christoffel  method  sketched  above  one  may  compute 
the  acoustic  and  piezoelectric  properties  for  wave  propagation  in  any  dir- 
ection. For  cubic  crystals,  three  directions  are  of  primary  importance 
from  the  point  of  geometry.  These  are  the  [100]  direction  (along  the  cube 
axis),  the  [no]  direction  (along  the  face  diagonal;  normal  to  the  dodeca- 
hedral planes  (110),  and  the  [m]  direction  (along  the  body  diagonal;  nor- 
mal to  the  octahedral  planes  (111)  ).  These  principal  directions  are  not 
necessarily  the  directions  of  greatest  interest  as  far  as  device  designs  are 
concerned,  but  they  do  bound  the  primitive  region,  the  irreducible  1 /48th 
of  the  unit  sphere,  for  these  crystal  classes.  The  region  is  shown  in  Fig- 
ure 1,  in  polar  coordinates,  in  terms  of  the  angles^  and  6 ; the  constant 
Q paths  are  portions  of  circles,  the  constants  paths  are  portions  of  radii, 
and  the  line  drawn  between  the  Iioo)  and  tin]  directions  has  the  equation 

a.yy  9 - sm  <j>  . (14) 

Numerical  values  for  Nm=ou/i,  |kj  and  | for  high  resistivity 
GaAs  are  given  in  Table  8 for  (100),  (110),  and  (ill)  cuts;  the  entries 
marked  by  dashes  indicate  eigenvector  degenaracies . For  wave  propagation  in 
the  [loo]  direction,  corresponding  to  the  (100)  cut,  the  modes  are  piezo- 
inactive. The  effective  elastic  constants  are  /c„,  £*•  , /c  • Cor" 

responding  to  the  (110)  cut  are  the  effective  elastic  values  c ^ 

^ Ol  ~ C ■£■  n + ^ I X.)  / 2.  j ^ nn  = ^ HH  * ^14  > A-*'1  d /Cc  = C it  tx.)  . 

The  piezoelectric  coupling  coefficients,  defined  by 

IU  = I ej  (15) 


are  zero  for  modes  a and  c;  for  mode  b the  value  is  IkJ  — I e,J/vXi^£  . 
As  seen  from  Table  8,  the  b mode  is  polarized  along  l°°  i]  , the  a mode  along 
[no]  and  the  c mode  is  perpendicular  to  both;  the  [no]  axis  has  particle 
motion  either  strictly  along  the  phase  progression  direction,  or  perpendi- 
cular to  it,  and  is  thus  an  "acoustic"  axis.27 
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TABLE  6.  ELASTIC  PROPERTIES  OF  SELECTED  II-VI  and  III-V  CRYSTALS  IN 


TABLE  7.  EXPERIMENTAL  VALUES  FOR  PHYSICAL  CONSTANTS  OF  GaAs  AT 


ROOM  TEMPERATURE. 


Quan- 

tity 

Unit 

Numeric 

Refer- 

ence 

- 

cn 

10+9Pa 

118.77 

16 

C 1 2 

53.72 

C44 

59.44 

-6 

Tcn 

10  /K 

-152. 

22 

Tc12 

r>. 

CO 

C\J 

1 

Tc44 

-37.8 

el4 

C/m2 

-0.16 

23 

£11 

pF/m 

116.7 

24 

en(T)/eo 

- 

12.73(l+pT) 

24 

(S=1.2x10"4/K 

P 

10+\g/m^ 

5.3169 

16 

1 


TABLE  8.  ELASTIC  AND  ELECTROMECHANICAL  PROPERTIES  OF  (100),  (110), 
AND  (HI)  CUTS  OF  GALLIUM  ARSENIDE. 


Quantity 

(100) 

(no) 

(111) 

i 

Cut 

Cut 

Cut 

2.363  2.617 
1.672  1.672 
1.672  1.237 


0 0 

0 6.08 

0 0 

1 1//T 

0 1//T 

0 0 

0 0 

0 
1 

0 1/  /2 

1/  JT 


0 


Wave  propagation  along  [nil  is  characterized  by  modal  elastic  stiff- 
nesses 

/C^=  ( /C„  + Z /Ca  + M + H /f)/3  , /Cb*/Cc  = C >C||  + ~ u ) • 

The  coupling  factor  for  mode  a is  |ka.l=  U &iJ//3£4  £ , while  the  other 
modes  have  zero  piezoelectric  coupling.  This  direction  is  also  an  "acoustic" 
axis. 


Cuts  of  GaAs  intermediate  between  the  three  principal  directions  are  # 
considered  next.  First,  the  great  circle  path  from  D100] ±0  [no]  ur>4h  e*o, &*+***> 
denoted  path  I,  is  described.  Figures  2,  3,  and  4 give  N»,  IkJ,  | */*,!, 
respectively,for  path  I.  The  degeneracy  in  the  b and  c modes  that  existsat 
Caoo]  is  broken  as  soon  as  this  specific  direction  is  departed  from.so,  for 
the  path  chosen,  the  ItLJ  may  be  extrapolated  back  to  the  end  point.  For 
this  path,  therefore,  one  has  libl-  (0,0,1)  ; this  vector  does  not  change 
during  path  I.  | | is  the  same  as  | with  a renumbering  of  the  com- 

ponents. 

Path  II consisting  of  the  great  circle  path  with  angles  4>  =45°, 

0°  ^ 9 $ sirT'O/VS)  , produces  the  results  shown  in  Figures  5,  6,  7,  and  8. 

This  path  starts  at  [110]  and  ends  at  [ml  ; the  latter  degeneracy  is  re- 

moved, as  described  above,  by  extrapolation.  For  mode  c,  the  components  do 
not  change:  | ^ 1 = ( ' /vT,  1 A/I , o)  . 

Path  III,  the  qreat  circle  path  from  [i 00]  to  [-i'll]  , has  angles  of 
0°  < 5*,  e-  sin^).  Figures  9,  10,  and  11  give  the  values  for  N„, 

IkJ  and  l^o.!  for  this  path.  The  behavior  of  for  path  111  is  exactly 

as  that  of  ) gfbl  for  path  II,  shown  in  Figure  8,  and  KJ  for  path  III  is 

the  same  asl'T*)  for  path  I when  the  1 and  2 components  are  interchanged. 

General  doubly  rotated  cuts  of  h!ah  resistivity  GaAs  are  considered  in 
Figures  12,  13,  14,  and  15.  Figure  12  gives  the  frequency  constants 
versus  Q for  various  <P  values.  One  sees  the  b and  c mode  degeneracy  at 
9 - r ' ( 1 / 1/3 ) , corresponding  to  [111.].  Of  interest  is  the  behavior  of  the 

c mode  between  2&. ia*  and  2e.110  , where  the  value  of  Ne  is  virtually  in- 

dependent of  ; this  appears  to  be  an  accidental  occurrence. 

Figures  13,  14,  and  15  show  the  piezoelectric  coupling  factors,  compu- 
ted from  Equation  (15),  for  modes  a,  b,  and  c,  respectively.  The  a mode 
has  its  maximum  value  at  <*>•*+  s',  e * s.rf’C , at  the  (111)  cut. 

The  b mode  value  is  a maximum  for  the  (110)  cut  at  this  drops 

off  as  9 is  increased,  and  becomes  zero  at  siVf'O/t/i')  and  beyond.  The  c mode 
behaves  in  a complementary  manner,  with  zero  coupling  for  <fi  = «s°  and 
& «.  , rising  beyond  this  9 value. 

Additional  details  relating  to  acoustic  wave  propagation  in  cubic  and/ 
or  other  crystal  classes  are  given  in  References  28  to  44,  and  69-78. 
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PIEZOELECTRIC  COUPLING  FACTORS  OF  DOUBLY  ROTATED  CUTS  OF  GALLIUM  ARSENIDE:  MODE 


MODE  b. 


VIBRATIONS  OF  DOUBLY  ROTATED  CUBIC  CRYSTAL  PLATES 


The  simple  thickness  vibrations  of  a crystal  plate  are  those  that  do  not 
depend  upon  lateral  coordinates.  Lateral  variations  make  the  problem  much 
more  compl icated1'5.  Taking  the  plate  thickness  direction  to  coincide  with 
the  wave  propagation  direction  discussed  in  the  previous  section,  the  solu- 
tion to  the  plate  vibration  problem  is  obtained  by  applying  the  appropriate 
boundary  conditions  at  the  plate  surfaces.  This  was  first  done  consistently 
by  Tiersten46* 47  for  traction-free  surfaces  and  an  insulating  crystal.  Be- 
cause of  the  piezoelectric  effect,  all  three  thickness  modes  are  coupled  at 
the  boundaries  in  the  general  case  of  complete  anisotropy.  Mass-loading  was 
considered  by  Yamada  and  Niizeki48,  who  derived  the  input  admittance  of  the 
plate  vibrator  as 


/ 3 i / 

X.puf  [i-2  UxM/xJ, 

(16) 

where 

C„  = e A Ah, 

07) 

= OJ  In  / , 

(18) 

and  kwis  obtained  from  Equation  (15).  C0  is  the  static  capacitance,  A is 
the  area  of  the  plate  portion  under  consideration,  and  2h  is  the  plate  thick- 
ness. Conditions  under  which  a single  mode  may  be  considered  by  itself  are 
discussed  in  Reference  6;  since  these  conditions  are  nearly  always  fulfilled 
for  the  materials  in  Table  6,  we  consider  a single  mode  in  the  sequel.  When 
the  conditions  are  not  fulfilled,  the  treatment  of  Reference  6 must  be  fol- 
1 owed . 

As  discussed  in  the  last  section,  if  the  crystal  is  not  insulating,  but 
semiconducting,  the  dielectric  permittivity  will  become  complex.  From  the 
fact  that  the  piezo-stiffening  term  then  becomes  complex,  becomes  complex 
as  well,  so,  from  Equations  (5),  (15),  (17)  and  (18),  the  various  factors  in 
Equation  (16)  for  the  input  admittance  of  the  crystal  vibrator  are  complex. 

Explicit  expressions  for  various  quantities  of  interest  are  given  in 
Table  9,  neglecting  carrier  diffusion,  where  Equation  (4)  can  be  written; 
the  effect  of  inclusion  of  the  diffusion  term  is  given  in  Reference  4. 

Viscous  loss  is  accounted  for  by  adding  to  the  elastic  stiffnesses  a 
term  , where  % is  a modal  viscosity.  This  has  an  effect  similar  to 

the  introduction  of  the  imaginary  permittivity  term,  but  the  frequency  vari- 
ation is  different.  One  distinguishes  between  the  acoustic  time  constant 

T„  = /X  , (19) 

and  the  electrical  time  constant 

V = £ /c r (2°) 
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TABLE  9.  PHYSICAL  PARAMETERS  FOR  INSULATING  AND  SEMICONDUCTING 
CRYSTALS. 


Quanti ty 


Insulator 


Semiconductor 


Symbol  Relation  Symbol  Relation 


by  the  different  frequency  dependencies  of  the  terms  containing  these  factors. 
For  example,  the  expression  for  the  effective  elastic  constant  in  the  presence 
of  both  acoustic  and  electrical  conductivity  losses  is 


= Z • [ i + i + k.V(coZ))]. 


When  the  plate  surfaces  are  coated  with  a thin  layer  having  a mass  m 
per  unit  area  on  each  surface,  then,  providing  that  the  coating  is  small 
enough  to  neglect  acoustic  transit  times  through  it,  the  effect  may  be  ex- 
pressed in  its  influence  on  Y-jnput  as  follows: 


C°/  { ^ * f (' 


lX  X. 


In  Equation  (22),  the  quantity  y is  the  reduced  mass  loading: 


yOL  = wi  / p h . 


The  relation  (22)  holds  also  when  the  terms  become  complex;  this  fact  will 
be  exploited  in  the  sequel.  An  air  gap  may  be  used  to  excite  the  resonators, 
rather  than  having  electrodes  deposited  directly  on  the  plates.  In  this 
case  the  effect  is  easily  taken  into  account  by  dividing  C0  and  the  in 
Equations  (16)  and  (22)  by  the  factor  ( -1  + c„  /cl)  , where  CL  is  the  air 
gap  capacitance  (or  value  of  a discrete  capacitor  placed  in  series  with  the 
crystal ) . 


EQUIVALENT  NETWORK  FOR  SEMICONDUCTING  PIEZOELECTRIC  PLATES 


From  the  input  admittance  given  in  the  last  section  the  equivalent 
electrical  network  can  be  synthesized.  This  will  be  briefly  sketched  here 
for  a single  mode,  although  the  network  can  be  extended  to  two  or  three  modesb. 
Starting  from  Equation  (16),  written  for  a single  mode,  the  admittance  of  a 
shunt  capacitance  C0  is  subtracted  out,  leaving 


i CJ  Ce  ( k*T)  / (l-k lT)  . 


Further  subtraction  of  a negative  series  C0  leave  the  admittance 


In  Equations  (24)  and  (25)  T equals  ‘tanx/x  . The  admittance  of  Equation 
(25)  is  realized  as  a transformer  shunting  a transmission  line,  the  trans- 
former representing  the  piezoelectric  transduction  mechanism,  and  the  line 
representing  the  mechanical  vibrations.  The  transformer  turns  ratio  is 


n = e. 


A /2h 


* 


and  the  transmission  line  admittance  is 

Y0  = 1 / A (*  ir»*.  . (27) 

When  the  conductance  is  not  zero,  the  network  appears  as  shown  in  Fig- 
ure 16.  Here  the  capacitances  CQ  are  shunted  by  conductances  G0,  where 

Ge  = <r  A / 2 h } (28) 

and  <r  is  tne  sample  conductivity  (reciprocal  resistivity).  When  diffusion 
effects  are  considered  one  will  have  an  effective,  frequency-dependent,  value 
of  <r  instead.  In  Figure  16,  Y0  and  will  also  be  complex  as  discussed 
above.  The  values  for  these  are  given  in  Table  10. 

Mass  loading  effects  are  taken  into  account  in  the  network  of  Figure 
16  by  simply  adding  in  series  with  the  transmission  line,  at  the  end  where 
the  piezoelectric  transformer  is  located,  an  inductance  L,  where 

L = * A . (29) 

This  relation  follows  from  performing  a network  synthesis  on  Equation  (22). 

The  effect  of  large  viscous  loss  on  the  behavior  of  the  equivalent  cir- 
cuit of  a vibrator  with  zero  conductivity  is  shown  in  Figures  17  and  18. 

In  Figure  17  is  shown  the  impedance  circle  representation;  2iM  = d-/Yinput. 
Figure  18  shows  the  real  and  imaginary  parts  of  the  input  impedance  as  func- 
tion of  frequency  normalized  to  the  frequency  where  the  lossless  transmis- 
sion line  would  be  one  quarter  wavelength  long. 

Measurement  methods  pertaining  to  low-Q  vibrators  are  discussed  in 
References  49-54.  With  the  inclusion  of  conductivity  terms  the  general  ap- 
pearances of  Figures  17  and  18  will  remain  as  shown,  but  the  expressions 
for  the  real  and  imaginary  parts  of  the  impedance  or  admittance  will  be  more 
complicated. 


NONDESTRUCTIVE  EVALUATION  OF  OXIDE  LAYER  MASS 


The  properties  and  performance  of  planar  integrated  circuit  devices55 
depend  on  the  production  of  inert  dielectric  oxide  layers  on  semiconductors. 

A gate  oxide  layer  may  typically  be  1000  J,  whereas  an  isolation  oxide  for 
passivation  is  perhaps  a micrometer  thick.  It  is  of  technological  interest 
to  have  a simple  and  rapid  nondestructive  test  of  the  amount  of  material 
added  to  the  semiconductor  in  both  cases.  For  this  purpose,  use  of  the  cry- 
stal as  an  acoustic  resonator  enables  the  critical  frequencies  associated 
with  the  resonance  to  be  measured56.  Since  these  may  be  measured  with  great 
accuracy  in  an  air  gap  holder57,  and  since  the  frequencies  are  simply  related 
to  the  parameters  of  the  equivalent  network,  the  mass  loading  due  to  the 
oxide  or  passivation  layer  may  be  deduced. 

In  the  case  of  gallium  arsenide,  some  of  the  various  oxide  species 
that  can  exist  are  enumerated  in  Table  11,  with  their  room  temperature 
densities.  From  the  table  one  can  see  that  the  densities  are  not  too  far 


FIGURE  16.  EXACT  EQUIVALENT  NETWORK  FOR  A SINGLE  DRIVEN  MODE  IN  A PIEZOELECTRIC,  SEMICONDUCTING  CRYSTAL  PLATE 
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TABLE  10.  ELECTRICAL  NETWORK  PARAMETERS  FOR  INSULATING  AND  SEMI- 
CONDUCTING CRYSTALS. 


Quanti ty 

Uni  t 

Insulator 

Semiconductor 

Symbol 

Relation 

Symbol 

Relation 

Static 

Capaci tance 

F 

Co 

6A/2h 

Co* 

Co+Go/ju 

Static 

Conductance 

S 

0A/2h 

Piezoelectric 

Transformer 

Ratio 

C/m 

n 

Ae/2h 

n* 

n 

Acoustic 

Admittance 

s/kg 

Yo 

1/A  pv 

V 

Y • F'k 

0 

Normal ized 

Frequency 

Variable 

_ 

X 

k h 

X* 

X • F* 

F = (l-k2/(l+jUT)} 


TABLE  11.  ROOM  TEMPERATURE  DENSITIES  OF  OXIDES. 


removed  from  that  of  GaAs,  averaging  about  3«io  kvA*3.  A knowledge  of  the 
particular  oxide  or  passivation  layer  consti tuent(s ) under  consideration  de- 
termines the  layer  thickness  once  the  added  mass  is  found  from  the  frequency 
shift. 

An  example  of  the  procedure  involved  follows:  The  sample,  in  plate 
form,  is  measured  to  determine  the  approximate  thickness.  With  the  sample 
orientation  known,  the  frequency  constants  are  known,  e.g.,  from  plots  like 
those  in  Figure  12.  The  corresponding  piezoelectric  coupling  factors  |k«,| 
are  likewise  known;  the  coupling  determines  how  strongly  the  mode  or  modes 
in  question  can  be  driven  by  the  measurement  apparatus.  The  frequencies 

"L  = (30) 

are  then  calculated;  about  these  frequencies  lie  the  resonance  regions.  Odd 
harmonics  of  are  also  piezoelectrically  excitable.  Once  the  resonance 
regions  are  approximately  determined,  the  sample  is  placed  in  an  air  gap 
holder,  and  the  gap  set  to  a prescribed  value.  Then  the  frequency  spectrum 
is  explored  by  any  of  several  methods58-63.  Some  point  on  the  immittance 
circle,  appropriate  to  both  the  measurement  instrument  and  the  resonator 
under  test,  is  selected  to  characterize  the  resonator  plate.  The  frequency 
at  this  point  is  recorded.  It  may  be  necessary  to  make  the  measurement  in  a 
modest  vacuum  for  best  results,  particularly  if  the  mode  a is  used.  The 
sample  is  then  removed  from  the  holder  and  processed  for  oxide  layer  growth 
or  passivation  layer  deposition,  etc.  After  completion  of  treatment,  it  is 
reinserted  into  the  air  gap  holder  and  remeasured  at  the  same  point  on  the 
immittance  circle  as  before  (for  example,  both  measurements  could  be  made 
at  the  condition  of  resonance,  where  the  imaginary  part  of  the  immittance 
is  zero  and  the  real  part  is  small;  or  it  could  be  made  at  the  frequency  of 
maximum  resistance,  or  of  minimum  impedance,  etc.). 

The  shift  of  the  chosen  critical  frequency  point  with  sample  treat- 
ment is  a very  sensitive  measure  of  p,  and  forms  the  basis  for  the  quartz 
micro-balance’'4-^".  In  the  present  i"r*2nce  the  frequency  shift  could  be 
correlated  with  p by  means  of  the  equivalent  circuit  of  Figure  16,  augmented 
as  described  by  the  addition  of  a lumped  inductor.  The  circuit  describes 
the  immittance  versus  frequency  of  the  vibrator,  and  the  effect  of  the  air 
gap  is  incorporated,  as  discussed,  by  modification  to  k^  and  CQ. 

Frequency  measurements  are  among  the  most  accurate  measurements  that 
can  be  made,  and  accuracies  of  fractions  of  a part  per  million  are  contnon- 
place.  To  get  an  appreciation  for  the  sensitivity  of  the  method,  consider 
a plate  having  a thickness  of  1 inn.  If  the  frequency  changes  by  lxl 0-6  on 
a normalized  basis,  and  if  this  frequency  change  was  due  solely  to  a change 
in  resonator  thickness,  then  the  thickness  change  would  be  1 mm  times  1x10"°, 
or  ioA  . Since  a fractional  frequency  deviation  of  1 x 1 0“ ^ ought  to  be  read- 
ily discernable  with  present-day  measurement  equipment,  and  since  the  plate 
thicknesses  are  apt  to  b'  only  fractions  of  a millimeter,  the  method  is  po- 
tentially capable  of  measuring  very  small  mass  changes.  The  method  does 
depend,  however,  on  the  accuracy  with  which  frequency  can  be  repeatably 
measured,  and  this  depends  on  a number  of  things.  Perhaps  chief  among  these  , 
aside  from  physical  reproducibility,  as  in  maintaining  the  air  gap  value,  is 


-5~r  — — — 


38 


temperature  variation.  For  quartz,  with  its  zero  temperature  coefficient 
cuts,  the  situation  is  eased.  Other  materials  lack  such  cuts,  and  this  in- 
cludes the  majority  of  crystals.  Approximate  values  for  the  temperature  co- 
efficients of  the  elastic  stiffnesses  of  gallium  arsenide,  calculated  from 
data  in  Reference  22,  are  given  in  Table  7.  Assuming  a (111)  cut,  the  first 
order  temperature  coefficient  of  frequency  for  the  a mode  is  approximately 
- 71  xio'Vk  ; the  value  for  a (110)  cut  (b  mode)  is  - 34  xicTVk  . One  sees 
from  these  numbers  that  a high  stability  oven  is  required  for  measurements 
in  the  10“®  range  for  gallium  arsenide  if  monolayer  masses  are  to  be  measured. 
A very  modest  oven,  however, # permi ts  one  to  measure  equivalent  thickness 
changes  in  the  100  to  1000  A range,  and  beyond.  This  is  just  the  regime 
encountered  in  integrated  circuit  electronics. 


CONCLUSIONS 


Acoustic  wave  propagation  in  cubic  crystals  has  been  considered  from 
the  standpoint  of  establishing  the  properties  germane  to  the  use  of  such  cry- 
stals as  piezoelectric  plate  vibrators.  An  equivalent  electrical  network 
was  developed  for  piezoelectric  semiconducting  plate  vibrators,  including 
the  presence  of  lumped  masses  on  the  plate  surfaces.  The  frequency  change 
due  to  the  added  surface  mass  was  shown  to  be  a very  sensitive  indicator  of 
the  quantity  of  mass  added,  and  can  serve  to  monitor  the  surface  layer.  Cal- 
culations were  carried  out  for  gallium  arsenide  doubly  rotated  plates;  these 
included  acoustic  velocities  and  piezoelectric  coupling  factors  for  the  three 
thickness  modes. 
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